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The ongoing experimental and theoretical effort aimed at understanding non-classical rotational 
inertia in solid helium, has sparked renewed interest in the supersolid phase of matter, its micro- 
scopic origin and character, and its experimental detection. The purpose of this colloquium is a) 
to provide a general theoretical framework for the phenomenon of supersolidity and 6) to review 
some of the experimental evidence for solid 4 He, and discuss its possible interpretation in terms 
of physical effects underlain by extended defects (such as dislocations). We provide quantitative 
support to our theoretical scenarios by means of first principle numerical simulations. We also 
discuss alternate avenues for the observation of the supersolid phase, not involving helium but 
rather assemblies of ultracold atoms. 
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The search for exotic phases of matter drives much of 
the fundamental research in condensed matter physics. 
Of particular interest are phases displaying simultane- 
ously different types of order. The archetypal exam- 
ple is the so-called supersolid, an intriguing, conjectured 
macroscopic manifestation of quantum mechanics, which 
has until now eluded direct, unambiguous experimental 
observation. In recent times, the search for this fasci- 
nating phase of matter has focused on solid helium, but 
alternate avenues to its possible observation are emerg- 
ing, notably in the area of cold atom physics. 

At the time of this writing, consensus is still lacking, 
as to whether the wealth of experimental evidence accu- 
mulated for solid helium (especially over the past eight 
years), indisputably points to the observation of the su- 
persolid phase of matter. Among the various proposed 
theoretical scenarios, none appears capable of accounting 
for all the puzzling, often conflicting experimental find- 
ings ( |Balibar[[2010 |. 

In this Colloquium, we make no attempt to offer the 
final word on this subject. Rather, we modestly limit 
ourselves to illustrating one theoretical interpretation of 
some of the existing helium phenomenology. Specifically, 
we shall describe a theoretical scenario in which extended 
defects, such as dislocations and grain boundaries, play 
a crucial role. Such a scenario is born out of, and sup- 
ported by, first principle microscopic calculations, based 
on Quantum Monte Carlo techniques. 

Motivated by the impetuous progress made in the 
context of ultracold atom physics, in light of the clear 
promise that such a field holds toward the realization 
of artificial quantum many-body systems, we discuss the 
possible observation of the supersolid phase in ultracold 
atom assemblies. In this case too, our theoretical pro- 
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posals are largely based on numerical evidence from com- 
puter simulations. 



A. Basic definitions 

What is a "supersolid" ? As it turns out, this con- 
cept defies in part one's intuition. Although a rigorous 
theoretical definition can be provided, it can be tricky at 
times to establish unambiguously whether a given physi- 
cal system meets all the requirements that warrant such 
a denomination. For this reason, we shall first proceed 
to state those conditions as clearly as we can, in order to 
set up the framework for the rest of our discussion. 

Taken literally, of course, the word supersolid means 
"superfluid solid". This implies that a supersolid sub- 
stance ought display the same basic qualities of a conven- 
tional solid, as well as those of a superfluid. This poses 
an intellectual challenge right away. For, one is hard 
put forming a mental picture of how the very different, 
and seemingly antithetical solid and superfluid physical 
behaviors, could be merged into a single homogeneous 
phase of matter, at least based on familiar macroscopic 
notions. Should one imagine a supersolid to be a peculiar 
solid, capable of flowing just like a liquid (and in addi- 
tion, without dissipation) ? Or, conversely, something 
more like a "stiff liquid" , resisting shear like a solid ? 
While there are some elements of truth in these tenta- 
tive descriptions, a rigorous, consistent definition of the 
supersolid phase is best achieved by taking on a micro- 
scopic approach, supplemented by the theoretical notion 
of long-range order. 

In order to keep the discussion as simple and trans- 
parent as possible, we shall for the moment restrict it 
to three-dimensional (3D) systems. In this case, both 
solid and superfluid phases have well-defined order pa- 
rameters that characterize them. In lower dimensional 
systems, solid and superfluid can be defined without ref- 
erence to order parameters (any such quantity vanishes 
identically at any finite temperature). Instead, one must 
speak of long-range (non-integrable) power-law correla- 
tions in the order parameter field (Kosterlitz and Thou 
less 1973). This difference, however, is not of fundamen- 



tal importance to our presentation, since modifications 
of theory required to deal with reduced dimensionality 
are relatively simple, and well understood. 



1. Order in a solid 

In a crystalline solid, elementary constituents (i.e., 
atoms or molecules) are spatially arranged in an orderly 
fashion, each occupying one of a discrete set of well- 
defined sites of a three-dimensional periodic lattice 1 . In 



order to make a more formal statement, let us introduce 
the (time-averaged 2 ) local density p(r) of particles form- 
ing the crystal. The average value of p(r) in a macro- 
scopic sample of crystal of volume Q is given by 
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d 3 r p(r). 



(1) 



Let 5p(r) = p(r) — p be the local deviation of the den- 
sity from the sample-averaged value p. In a system that 
does not break translation invariance, such as a gas or a 
liquid, p(r) — p and Sp(r) = 0. In a crystal, on the other 
hand, translational invariance is broken, i.e., 5p(r) does 
not vanish identically, and ordering is expressed through 
the following condition: 



6p(r) = 8p(r + T), 



(2) 



holding for any vector T belonging to the discrete set 
of lattice vectors (i.e., vectors connecting any two lat- 
tice sites) defining the crystal. Consider the three- 
dimensional Fourier transform of Sp(r) 



p(k) = 



d 3 r -- 4k " 



5p(r) 



(3) 



Condition ^ implies that p(k), and its squared magni- 
tude 5*(k), known as the Static Structure Factor, display 
peaks in correspondence of wave vectors G in the so- 
called reciprocal lattice, i.e., fulfilling the condition 



G • T = 2irn 



(4) 



for any vector T defined as above, with n being an inte- 
ger. These peaks in S(k) manifest themselves experimen- 
tally as maxima in the intensity of light scattered off a 
crystalline sample, at specific angles ( Chaikin and Luben- 
sky 1995 ) . Because condition ^ is assumed to hold for 
a macroscopic sample, one speaks of density long-range 
order (LRO). 

Broken translation invariance has the immediate 
macroscopic consequence of imparting to the system the 
well-known property of resistance to shear, which a fluid 
does not possess. It is important to note that, while LRO 
implies broken translation invariance, the converse is not 
true. Glasses, for example, are substances which share 
many of the macroscopic properties of crystalline solids, 
in which translation invariance is broken, but Sp(r) fluc- 
tuates disorderly and no LRO exists. As we shall discuss 
later in this chapter, there is in principle nothing prevent- 
ing some glassy systems from turning superfluid; how- 
ever, the denomination supersolid does not apply here, 
due to the lack of density LRO. In these cases, one ought 
more properly speak of superglass. 



1 We assume for simplicity a Bravais lattice. 



2 In any real crystal, thermal and quantum fluctuations cause the 
local density to deviate instantaneously from its time-averaged 
value. 
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Another point that need be made at the outset, is that 
condition ^ occurs in a solid spontaneously, exclusively 
as a result of interactions among elementary constituents, 
at specific thermodynamic conditions (i.e., temperature 
and pressure). Such a spontaneous breaking of trans- 
lation symmetry is an integral part of the definition of 
supersolid. One could certainly imagine externally im- 
posing a density modulation to a superfluid liquid. For 
instance, this could be accomplished by adsorbing few 
layers of 4 He on a strongly at t ractiv e substrate, such as 
graphite (Crowell and Reppy 1993). The local density 
of successively adsorbed 4 He layers inevitably reflects the 
corrugation (i.e., LRO) of the underlying graphite sub- 
strate, to lesser a degree for layers increasingly removed 
from it. Referring to any superfluid adsorbed 4 He layer 
as "supersolid" , on account of the density modulation 
imposed by a corrugated substrate, is fundamentally in- 
correct and constitutes a misnomer. 



2. Order in a superfluid 

Superfluidity is the property of a substance sustain- 
ing persistent, dissipation-less flow. It was first observed 
in the liquid phase of the most abundan t isotope of he- 



lium ( 4 He), independently by Kapitsa (Kapitza 1938 ) 



and Allen and Meisener ( Allen and Misener[ 1938) in 
1938. Its fundamental relevance extends considerably be- 
yond the physics of condensed helium. For example, it is 
now widely accepted that superconductors are essentially 



charged superfluids (Tilley and Tilley 1990). 



A useful phenomenological model introduced by Tisza 
shortly after the discovery of superfluidity in liquid 4 He 
( Tisza 1938 ) , formally expresses the local density p(r) of 
a superfluid system, as a sum of two contributions: 



p(r) = p s (r) + p N {r) 



(5) 



where the subscripts S and N stand for "superfluid" and 
"normal" component. The superfluid component carries 
no entropy and can sustain flow indefinitely, whereas the 
normal component is subjected to dissipation, as any reg- 
ular fluid. Henceforth, we shall use the notation ps and 
Pn to refer to values averaged over the whole system. 
In a translation invariant system, it is pjv(r) = pn and 
ps(r) = ps. Superfluid behavior typically sets in through 
a continuous phase transition at a temperature T c , below 
which ps becomes nonzero, and increases monotonically 
as the temperature T — > 0. In a translation invariant 
system, ps{T = 0) = p, i.e., the average density of the 
system. Phrased differently, the superfluid fraction ps / p 
approaches unity as T — > 0. On the other hand, Leggett 
showed (Leggett 1970) that, for any superfluid system 



breaking translation invariance, it is ps(T = 0) < p. 

A wealth of theoretical and experimental work, span- 
ning now over seven decades, has afforded satisfactory 
theoretical understanding of the microscopic origin of su- 
perfluidity, in 4 He and other systems. In particular, it is 
now understood that superfluidity is a macroscopic man- 



ifestation of quantum particles behaving collectively as 
classical complex fields (iKosterlitz and Thoulessl 11973 



Svistunov 19921, and that it is crucially underlain by 



a type of quantum statistics (Bose) since description in 
terms of classical fields emerges in the limit of large oc- 



cupation numbers (Glauber 1963 Langer 1968 1969). 



In three dimensional systems, superfluidity is accom- 
panied by Bose-Einstein condensation (BEC) ( |Leggett| 
2006), a collective phenomenon that takes place at low 
temperature. Indeed, there is virtually unanimous agree- 
ment that, if superfluidity is to occur in a system of 
identical particles obeying Fermi statistics, some physical 
mechanism must be present, allowing for the formation 
of pairs (e.g., electronic Cooper pairs in superco nductors, 
or atomic pairs in 3 He (Tilley and Tilley 1990) and cold 
atomic vapors). These composite objects act in some 
sense as bosons, and as such can undergo BEC 3 . 

BEC consists of the occupation of just one quantum- 
mechanical single-particle state by a finite fraction of all 
N particles in the system 4 . In order to be regarded as 
such, BEC must take place in a macroscopic sample of 
matter; mathematically, this is expressed by taking the 
thermodynamic limit, i.e., N, f2 — > oo, but the density 
p = N/n remains finite. 

In a many-particle system that enjoys translation 
invariance, like a liquid, the single-particle quantum- 
mechanical state into which particles can condense, is 
that of a free particle of momentum hk. = 0. In order to 
render the above statement more quantitative, we intro- 
duce the momentum distribution 



n(k) = (^(k) i>(k))/N 



(6) 



where (...) stands for thermal expectation value and ^>(k), 
■0^(k) are Bose annihilation and creation operators of a 
particle of momentum Kk. The operator %jv (k)^(k) is 
the number operator for particles of said momentum. 
In a Bose condensed system, n(k) will take the form 



n(k) = n <5(k) + n NC (k) 



(J) 



The first term refers to the condensate, and the quan- 
tity n is referred to as condensate fraction. It should 
not be confused with the superfluid fraction, namely the 
fraction of the system that can flow without dissipation. 
These are two conceptually distinct quantities. In par- 
ticular, n a is the square of the amplitude of the complex 
order parameter for the superfluid. The second term of 
([7]) represents the contribution to the momentum distri- 
bution coming from states of nonzero momentum, none 



3 The relationship between superfluidity and BEC is straightfor- 
ward only in 3D; formally BEC is not required for superflu- 
idity, as evidenced by finite temperature superfluidity in two- 
dimensional (2D) systems, where BEC is replaced by topological 
order and power-law correlations. 

4 For simplicity, but with no loss of generality, we restrict our 
discussion to systems of bosons of spin zero. 
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of which is occupied macro scopically (i.e., by a finite frac- 
tion of all the particles in the system). The condensate 
fraction approaches a value close to 100% at low temper- 
ature, in a weakly interacting Bose gas. Interactions have 
the effect of depleting the condensate; for example, the 
most current experimental estimate of n in liquid He 
at zero temperature and at saturated vapor pressure, is 
~ |2000l ). 



about 7.5% (iGlyde et al. 



In order to elucidate how a finite condensate frac- 
tion corresponds to a kind of order, consider the Fourier 
transform of the momentum distribution, known as one- 
particle density matrix 



n(r,r') = (^(r)ip(r')) 



(8) 



Here, the field operators ^(v) and ip ( r ) annihilate and 
create a particle at position r. They are Fourier trans- 
forms of ip(k) and V^(k). 

In a system enjoying translation invariance, such as a 
liquid or a gas, it is n(r,r') = n(r — r'). If translation 
invariance is broken, as in a crystal, it is customary to 
consider, for simplicity, the spatially averaged function 

n(r) = ^J d 3 r'n(r',r' + r) . (9) 

It can be easily shown that condition ^ is equivalent to 

n(r) — > n a as r — > oo , (10) 

implying that a state of a system comprising a thermody- 
namic number of particles, in which a particle is removed 
at a given position r, has a finite quantum-mechanical 
amplitude over a state in which an identical particle is 
removed at an arbitrarily large distance away from r. 
Condition ( 10 1 expresses the presence of a type of order, 



referred to as off-diagonal long range order (ODLRO). 

One way to come to terms with the mind-boggling no- 
tion of ODLRO, consists of accepting that, because par- 
ticles are indistinguishable and can trade place with one 
another, each particle can be seen as being essentially 



delocalized throughout the whole system (Kohn 1964 
Leggettl |1970h. 



3. Order in a Supersolid 

Henceforth, we shall adopt the term "supersolid", to 
refer to a homogeneous 5 phase of matter in which both 
LRO and ODLRO, as defined above, exist simultaneously 
and appear spontaneously for the same species of parti- 
cles. The definition of ODLRO is the same in a system 
featuring crystalline order. Conventional momentum is 
replaced by crystal momentum, but BEC, if it occurs, 



still takes place in the k = state (up to a reciprocal 
lattice vector), by symmetry. Moreover, since a crys- 
tal breaks translation invariance the superfluid fraction 
ought to saturate to a value less than unity as T — » 0. 

At first, one might think that the intrinsic localiza- 
tion of particles characterizing a crystal phase, and the 
derealization ensuing from the condition ( 10 1 (as well 



as superfluidity) ought to be incompatible. Indeed, this 
was the conclusion reached by Penrose and Onsager, ar- 
guably the first to pose the question of the existence of 
a supersolid in 1956. In their seminal paper, they con- 
tended that no ODLRO can exist in a crystalline solid. 
Their argument is based on a variational model of a 
crystal with atoms "pinned" at equilibrium lattice po- 
sitions (Penrose and Onsager 1956). The statement is 



that particle localization, due to crystallization or other 
causes (e.g., disorder), acts so as to prevent macroscopic 
quantum-mechanical exchanges of indistinguishable par- 
ticles, thereby de facto removing all effects of quantum 
statistics (including BEC). 

While seemingly plausible, and in fact re-iterated even 
many years later ( |Nozieres 1995 1, the argument of Pen- 
rose and Onsager is no formal proof of the non-existence 
of a supersolid. Indeed, shortly after the publication of 



the article by Penrose and Onsager, Gross ( Gross 1957 



1958 



) showed that a superfluid system described by the 
non- inear classical field equation may feature a density 
wave modulation. This finding, which qualifies as the 
first theory of the supersolid phase, was essentially over- 
looked, in part due to incomplete understanding of the 
conditions leading to the validity of the calculation by 
Gross. In 1962, C. N. Yang (Yang 1962) proposed that 



in a crystal of helium, in which atoms enjoy a high de- 
gree of quantum derealization, exchanges of adjacent 
atoms might in fact be significant. Successively, Leggett 
(Leggett 1970) revisited Yang's argument, suggesting an 



experiment in which the (presumably small) superfluid 
response of a 4 He crystal could be measured. 

A different, intriguing scenario for the occurrence of 
supersolid order in a highly quantum crystal, came from 
Andreev and Lifshitz (Andreev and Lifshitz 1969]), and 



Chester (Chester 1970), and is based on the predicted 



As opposed to simple coexistence of two phases, each featuring 
only one kind of order. 



high mobility, at low temperature, of point defects like 
vacancies or interstitials, which are thermally activated 
in most materials, but may be present even in the ground 
state of such quantum many-body systems as 4 He. These 
delocalized defects are often referred to as "zero-point 
defects" , to emphasize the notion that they are not pro- 
duced by deliberately removing (adding) atoms from (to) 
the lattice. Instead, the lattice period and the unit cell 
volume, Oo, are automatically adjusted to have a non- 
integer number of particles per unit cell, pfin, even in the 
fully equilibrated state. Conceivably, a gas of such repul- 
sively interacting point defects ought undergo BEC, and 
turn superfluid, at low temperature. More recently, it 
was pointed out that the "zero-point defects" supersolid 
scenario is the only one possible in perfect continuous 
space crystals (Prokof'ev and Svistunov 2005) because 
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superfluid states are necessarily gapless with respect to 
adding and removing particles to/from the system. Con- 
sequently, integer pQ in a supersolid may occur only as 
a zero-measure coincidence. 

A crystal of 4 He has always been regarded as the most 
likely candidate to display supersolidity, arising through 
any of the above scenarios, even though no reliable quan- 
titative predictions have ever been furnished of the mag- 
nitude of the effect, or the transition temperature. This 
is perhaps the main reason why, despite an intense exper- 
imental effort, no unambiguous experimental observation 
of supersolid behavior in 4 He was reported for decades, 
following the original theoretical predictions. In 1992, a 
rather sobering assessment of the prospects of successful 
observation of a supersolid phase of helium, was offered 



by Meisel (Meisel 1992) 



An exciting, unexpected turn of events took place in 
2004, with the discovery, by Kim and Chan, of a down- 
ward shift in the period of a torsional oscillator filled 
with solid 4 He, at temperatures below 250 mK (Kim 



and Chan 2004a|b ). The interpretation of these find- 



ings in terms of a superfluid mass decoupling within the 
helium solid unquestionably appears legitimate, serious 
and plausible. However, as we shall discuss in the next 
sections, the final word is yet to be spoken as to whether 
the effect first reported by Kim and Chan, and since then 
confirmed by a number of other groups, indeed signals the 
onset of a new phase of matter. 

Irrespective of the above, though, there is no question 
that the discovery has sparked a renewed interest in the 
phenomenon of supersolidity, from which novel, interest- 
ing lines of research have spun. For one thing, achieving 
a complete, thorough understanding of the low tempera- 
ture phase diagram of solid helium has proven a surpris- 
ingly complicated, and surely intriguing problem on its 
own, both from the theoretical and experimental stand- 
points. Many puzzling, and often apparently conflicting 
experimental findings, have cast doubts on the applica- 
bility to solid 4 He of microscopic pictures such as that of 
Andreev, Lifshitz and Chester. In particular, significant 
attention is being directed to inhomogeneous scenarios, 
involving such extended defects as grain boundaries, as 
well as dislocations. 

Secondly, the search for the supersolid phase of mat- 
ter has now extended beyond solid helium, in a new 
and exciting direction. The scientific breakthroughs in 
cold atom physics that have characterized the past two 
decades, make it now feasible to investigate novel phases 
of matters, not only by providing remarkably clean and 
controlled experimental many-particle systems, but also 
by allowing one to fashion artificial inter-particle poten- 
tials, not arising in any known condensed matter system. 
The aim here is twofold: on the one hand, the high degree 
of control afforded by cold atoms, may render the unam- 
biguous identification of the sought phase more direct, 
and easier than in solid helium. More fundamentally, the 
possibility of "tweaking" the inter-particle potential per- 
mits one to address basic questions, such as: which kind 



of interaction, if any, can underlie a supersolid phase 



B. Experimental detection 

Because it is a manifestation of quantum mechanics 
on a macroscopic scale, one would be inclined to as- 
sume that the observation of supersolid behavior should 
be relatively straightforward and unequivocal, the only 
issue being whether the relevant phenomenology occurs 
at experimental conditions (e.g., temperature) that are 
presently accessible. Observation of superfluidity in liq- 
uid He is made relatively simple by the stunning effects 
that can be produced in laboratory, all of which exploit 
the unusual flow of the superfluid liquid. The situation 
is more complicated in the case of a supersolid, which is 
expected to retain the mechanical properties of a regular 
solid, notably resistance to shear, preventing it to flow 
like a liquid 6 . 

Strictly speaking, some demonstrations of superfluid 
behavior, such as persistent flow, or the well-known foun- 
tain effect (Tilley and Tilley 1990), should be feasible in 



a solid as well, with the important caveat that a super- 
solid will support dissipation-free flow only of its own 
particles; with respect to any other external object, it 
will behave as a conventional solid. Clearly, however, 
the ease and spectacular immediacy of the observation 
of such effects, when they take place in the liquid super- 
fluid phase, may be difficult to replicate in the solid. In 
general, therefore, the supersolid nature of a condensed 
matter system can only be typically ascertained by means 
of relatively indirect observations. 

Another potential problem lies in the fact that, as men- 
tioned above, in the low temperature limit neither the 
superfluid nor the condensate fractions of a supersolid 
will approach 100%, due to interactions among particles 
and the lack of translation invariance. Rather, they will 
both saturate to values that could conceivably be very 
small, rendering the measurement of any "superfluid sig- 
nal" problematic. 



1. Non-classical Rotational Inertia 

Perhaps the simplest way of detecting experimentally 
supersolid beh avior, was suggested by Leggett in 1970 
(Leggett 1970 1, and consists of measuring a reduction of 



the moment of inertia of a solid sample, with respect to 
its classical value. This effect is commonly referred to 
as Non-Classical Rotational Inertia, or NCRI. Consider 
a cylindrical vessel (cell), filled with a given substance, 
and let the cell be connected to a torsion rod, allowing 



6 Flow arising as a result of plastic deformation of a supersolid, 
is expected to be radically different from that in a regular solid. 
The corresponding theory, however, is still far from being com- 
plete (see also below) 



it to execute small harmonic oscillations about its axis, 
as shown in Fig. [TJ This experimental set-up is typically 
referred to as Torsional Oscillator (TO). 

The measurable resonant period of oscillations is pro- 
portional to the square root of the moment of inertia of 
the system, which of course includes a contribution from 
the vessel itself, as well as from the substance enclosed 
in it. 



Torsion rod 



Torsion cell 
containing 
helium 

Drive - 




Detection 



FIG. 1 Schematic of the experimental Torsional Oscillator 
set-up aimed at measuring non-classical rotational inertia. 
The cell is filled with a sample of condensed matter, whose 
superfluid fraction affects the moment of inertia of the overall 
system. 




FIG. 2 Color online. Schematic of an experiment measur- 
ing mass flow through a supersolid. Superfluid liquid flows 
across the central section filled with supersolid, of the same 
substance. The setup as shown can work only at the melting 
pressure. In the case of 4 He, modifications involving Vycor 
segments and temperature gradie nts allow one to extend ed 
measurements to higher pressures ( Ray and Hallock 2008 1 . 



For these reasons, careful, concurrent control measure- 
ments on similar, non-superfluid systems 7 are required, 
in order to ascertain by comparison whether what is be- 
ing observed is indeed a genuine superfluid transition. 

There are a number of puzzling observations, con- 
nected with the onset of NCRI, which have so far eluded 
a satisfactory explanation. We come back to this point 
below, when discussing experimental evidence of superso- 
lidity of 4 He that is available at this time, NCRI measure- 
ments constituting the bulk thereof. We shall propose a 
possible theoretical interpretation of (at least some of) 
the phenomenology in that context, in Sec. In] 



If below a given temperature, the resonant period 
shows a monotonic, continuous decrease, one way to in- 
terpret this is to assume superfluid mass decoupling from 
the rotation of part of the sample inside the vessel. In 
other words, a fraction thereof (normal) keeps oscillating 
together with the apparatus, while the remaining (super- 
fluid) part remains at rest, therefore not contributing to 
the rotational moment of inertia. Formally, 



I(T) = I class (l-p s (T)/p) 



(11) 



where I c iass is the classical moment of inertia, regarding 
the whole sample as normal. Thus, a measurement of 
the resonant period is a measurement of the superfluid 
fraction of the system, as long as /class and the elastic 
properties of the entire apparatus remain constant as the 
temperature is lowered. 

The NCRI technique is particularly well-suited to 
probe possible superfluidity of a solid, as it does not 
require any extra circuitry to produce relative flow be- 
tween the normal and superfluid components. However, 
the NCRI-TO experiment is not without some ambiguity, 
in part precisely because it is indirect, i.e., it does not in- 
volve any actual measurement of mass flow out of the sys- 
tem. Indeed, a number of different physical mechanisms 
have been proposed, not involving a superfluid transi- 
tion of the helium crystal, that could conceivably cause a 
change in the oscillation period, thereby giving rise to a 
spurious NCRI signal ( Clark et al. 2006 1 . An example of 



2. Mass transport through a supersolid 

A different way of detecting supersolid behavior is 
based on flow, not of the supersolid, but through the su- 
persolid. Consider the experimental set-up schematically 
shown in Fig. [2] One wishes to establish flow of a fluid 
through a pipe, a section of which is filled with solid of 
the same substance. The most obvious example, one with 
which we are all familiar, is that of a water duct. During 
a cold Winter month, water can freeze in a section of the 
duct. As we know from experience, water flow (for ex- 
ample, in the direction shown by the arrows in Fig. |2j), 
ceases as soon as ice (solid water) clogs the duct. There 
is, of course, nothing peculiar about water in this regard; 
the same observation will be made if any other substance 
fills the pipe, as long as freezing takes place as shown. In- 
deed, not even superfluid liquid helium may flow through 
the same section of the pipe, if filled with normal solid 
helium. 

On the other hand, some flow will take place through 
supersolid helium. In other words, a superfluid liquid will 
flow through a supersolid made of the same elementary 
constituents. The above experiment allows one to gain 
intuitive, vivid understanding of the somewhat cryptic 



such a mechanism is given by a possible structural trans- 
formation taking place in the crystal at low temperature. 



7 For example, in the case of solid 4 He an obvious control experi- 
ment consists of performing the same measurement on solid 3 He, 
a Fermi system not expected to undergo a supersolid transition 
— certainly not at the same temperature as 4 He. 
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idea of ODLRO, as expressed by condition (10); for, one 
may think of helium atoms in the superfluid phase, flow- 
ing toward the frozen blockage (left part of Fig. [2]), es- 
sentially "re-emerging" on the right side of it, through a 
process of long exchanges with (indistinguishable) atoms 
in the supersolid. Thus, "flow" through a supersolid nec- 
essarily must involve identical particles, as it is underlain 
by quantum exchanges. An impurity, namely a particle 
distinguishable from those that make the supersolid, will 
not move through it any more than through a regular 
crystal. In this respect, a supersolid acts very differently 
from a superfluid, which allows a foreign particle to move 
through it with very little drag (no drag at all at T = 0). 

Experiments such as the one described above, have 
been performed on solid helium only relatively recently 
( [Ray and Hallock[ [2008| . They are also not without 



drawbacks, e.g., the difficulty of inferring a reliable value 
for the superfluid fraction. We shall describe the main 
results of these experiments, and the progress that they 
have afforded on the understanding of supersolidity in 
4 He, in Sec. |n| 



nal, make it all but out of question to address the issue 
of supersolidity of 4 He at the present time. Moreover, 
there exist more robust techniques (e.g., x-ray diffrac- 
tion) that allow one to establish unambiguously whether 
a condensed matter system of interest has crystal order; 
it therefore makes sense to carry out separate measure- 
ments, aimed at probing individually superfluidity and 
density LRO. Because in the presence of interactions, the 
superfluid fraction is usually greater in magnitude than 
the condensate fraction, there is generally a preference 
among experimentalists for measuring the former. 

In a different physical contexts, however, the measure- 
ment of n(k) is inarguably the most efficient way of es- 
tablishing the supersolid nature of an underlying many- 
particle systems. This is typically the case for assemblies 
of cold atoms, whose momentum distribution lends itself 
to direct imaging, by means of time-of-flight measure- 



ments ( Greiner et al. 2002 ) 



C. Theoretical studies 



3. Momentum Distribution 

The momentum distribution n(k) (Eq. [6]), can in prin- 
ciple expose in full the supersolid character of a system. 
For a crystal of "Boltzmannons" , namely distinguishable 
quantum particles, n(k) is a Gaussian, whose semi-width 
is determined by zero-point motion. Particle indistin- 
guishability and quantum statistics cause measurable de- 
viations of the momentum distribution from the Gaus- 
sian behavior. Specifically, the one-body density matrix 
n(r), which is the Fourier transform of n(k), decays expo- 
nentially in a normal (i.e., non-supersolid) Bose crystal, 
conferring to n(k) a much slower decay with k than the 
Gaussian one ( Boninsegni , 2009 ) . 

In the presence of ODLRO, n(r) settles asymptoti- 
cally on the finite value n at large r. In the case of 
a supersolid, however, it must also take on an oscillatory 
behavior, reflecting lattice periodicity. This is because, 
in a solid, long many-particle exchanges which underlie 
ODLRO involve particles situated preferentially at lattice 
sites. That means that n(k) in a supersolid will feature 
satellite peaks at reciprocal lattice vectors. 

Thus, in principle a measurement of n(k) could pro- 
vide the most direct assessment of whether a system is 
supersolid or not, as both types of order are probed simul- 
taneously. In practice, however, the feasibility of such a 
determination depends on the strength of the signal, i.e., 
7io i as well as on the inherent sensitivity of the measuring 
apparatus. 

For example, the momentum distribution of solid he- 
lium can be studied experimentally by deep inelastic neu- 
tron scattering. The resolution of the most advanced ex- 
perimental facilities available nowadays, as well as the 
uncertainty associated with the analysis of the raw data 
and the expected, sheer feebleness of the supersolid sig- 



From the theoretical standpoint, the investigation has 
proceeded along several fronts. This Colloquium is not 
aimed at providing an overview and/or critical compari- 
son of all the different theories. For example, we shall not 
review nor discuss any of the scenarios which attribute 
the observed NCRI to effects that are mechanical in ori- 



response of solid 



Yoo and Dorsey 



4 He ( 


Pratt et al. 


2009 


. Likewise, 



2011 



Reppy 



2010 



we shall not review 



the "vortex liquid" theory 8 (Anderson 



2007 



2008). All 



of these theoretical proposals are the subject of intense 
current debate, and interested readers should check the 
pertinent references. Instead, we shall restrict here our 
discussion to scenarios that attempt to connect experi- 
mental observations to first principle calculations, carried 
out on increasingly realistic models of a helium crystal. 



1. First principle calculations 

The expression "first principle" , in this Colloquium, 
means that one starts from a quantum-mechanical many- 
body Hamiltonian H for a system of N identical particles, 
namely: 



^ = ^^E V ' + ^ r i' r 2,...r J v) 



(12) 



8 It is worth mentioning that, while the "vortex liquid" theory, 
which assumes a superfluid ground state for an ideal 4 He crystal, 
is in direct contradiction with existing numerical evidence, it 
contains some of the same mesoscopic phenomenology present 
in the so-called Shevchenko state (Shcvchcnko 1987), discussed 
below. 



where m is the mass of each particle in the system, and 
where V is the potential energy of interaction. In all 
cases of interest here, V can be expressed as the sum of 
terms involving pairs of particles, i.e., 



V(ri,r 2 ,...rjv) =^«(i 



(13) 



'<3 



where the function u(rj, Tj) is known either from accurate 
ab initio quantum chemistry calculations (e.g., the inter- 
action between two helium atoms by Aziz et at (1979])), 
or, in the case of cold atoms, can be as simple and funda- 
mental as the potential of interaction between two elec- 
tric dipoles. 

The calculation of cogent physical quantities, char- 
acterizing equilibrium thermodynamic phases of many- 
body systems modeled as in Eq. [12J can be carried out 
to a remarkable degree of accuracy by means of numer- 
ical simulations, based on Quantum Monte Carlo tech- 
niques (see, for instance, Ref. |Boninsegni et al. ( 2006b l). 
This computational methodology is fairly mature, and 
thoroughly described in the literature; we therefore re- 
fer interested readers to the original references. What 
is important to stress here, is that these methods allow 
one to obtain direct, unbiased numerical estimates of all 
the quantities that are relevant in the study of the su- 
persolid phase of matter, namely the one-body density 
matrix n(r), the supcrfluid fraction ps(T) and the static 
structure factor (standard finite-size scaling procedures 
allow one to extrapolate results to the thermodynamic 
limit). Essentially all of the theoretical results presented 
and discussed in the following sections, are obtained by 
means of Quantum Monte Carlo simulations. 



2. Lattice Supersolids 

A different approach, also affording fundamental in- 
sight into the nature of the supersolid phase, the tran- 
sition to a supersolid from a normal solid or from a 
supcrfluid, and the role of point defects, makes use of 
much simpler models than Eq. |12[ wherein particles 
are confined to moving on discrete lattices. In this case, 
spontaneous breaking of translation invariance is defined 
with reference to the discrete translation symmetry of 
the Hamiltonian. The simplicity of the models utilized, 
and the possibility of carrying out very accurate numer- 
ical calculations, make it possible to address a number 
of theoretical issues relevant to the problem of super- 
solidity, in some cases offering a valuable guide to the 
experimental investigation. Indeed, some of the earliest 
microscopic theoretical investigations aimed at assessing 
the possible existence of a su persolid phase of 4 He, made 
use of lattice Hamiltonians ( |Liu and Fisher 1973 



Mat- 



suda and Tsuneto 1970). Moreover, while for a long time 
lattice models were regarded as little more than a useful 
theoretical device, impressive advances in optical lattice 
technology render it now possible to synthesize and in- 
vestigate in the laboratory artificial many-body systems, 



very accurately realizing the physics embedded in those 
models. Therefore, one may conceivably be able to carry 
out a direct comparison of theory and experiment, enjoy- 
ing a degree of accuracy not attainable in solid helium 
experiments. 

There exist, however, fundamental differences be- 
tween continuous space and lattice supersolids, which in- 
evitably limit the understanding that any one of them 
can afford of the physics of a system in continuous space. 
A lattice "crystal" is one in which basis vectors are linear 
combinations of those of the underlying lattice. Such a 
solid period cannot be varied continuously, i.e., a lattice 
crystal is by definition incompressible. Consequently, su- 
persolid phases can be obtained in a rather trivial way, 
simply by doping the insulating crystal with holes or par- 
ticles (provided that no phase separation occurs of the 
dop ed system into a superf l uid an d an undoped crys- 
tal (Batrouni and Scalettar 2000)). While continuous 



space crystals can lower their energy by eliminating va- 
cancies and interstitials, using complete atomic layers 
and adjusting the lattice constant accordingly, an anal- 
ogous mechanism is forbidden if crystalline sites must 
remain pinned to an underlying space lattice. Conse- 
quently, in lattice crystals point defects can be intro- 
duced "by hand"; their presence is not necessarily of 
"zero-point" origin, and their number is conserved. The 
second crucial difference between continuous and pinned 
lattice solids, is the absence in the latter of long-range 
elastic forces between point defects. Thus far, all lattice 
supersolids studied in the literature (in fact, we are not 
aware of a single counter-example) belong to the cate- 
gory of doped insulators, which is fundamentally different 
from the zero-point defect picture relevant for continuous 
space systems. Heretofore, we restrict our discussion to 
continuous-space supersolids, making use of lattice sys- 
tems occasionally, for illustration purposes only (e.g., in 
Sec. II. B). Readers interested in lattice supersolids are 



referred to the vast literature on this subject, of which 



Batrouni and Scalettar 


(200C 


I ; Boninsegni 


(2003) 


; Bonin- 


segni and Prokof'ev ( 


2005); Chen et al. 


(2008 


); Dang 


et al. 


(2008); Hebert et al. ( 


2001); Hcidan 


an and Damle 


(2005 


; Melko et al. (|2005l; 


Sengupta et aL|(|2005l; Wes- 


sel and Troyer (2005|) are just a representative sample. 



In the following two sections, we shall a) review the 
present status of the theoretical and experimental re- 
search on supersolidity in 4 Hc, and 6) discuss theoreti- 
cal predictions for supersolid behavior in systems of cold 
atoms. 



II. SUPERSOLIDITY IN 4 HE 

The experimental search for the supersolid phase of 
matter, has focused over the last four decades almost ex- 
clusively on the most quantum solid in Nature, whose el- 
ementary constituents obey Bose statistics, namely 4 He. 
Helium, unlike any other naturally occurring substance, 
remains a liquid under the pressure of its own vapor, all 
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the way down to zero temperature. It solidifies under 
some moderate pressure, (approximately 25 atm), form- 
ing a hexagonal close-packed (hep) crystal. It is worth 
noticing upfront two aspects, that complicate the unam- 
biguous detection of supersolid behavior in a helium crys- 
tal. The first is the resilience of superfluidity in the liquid 
phase, which is predicted to remain supcrfluid in the low 
temperature limit at arbitrarily high density, i.e., in the 



metastable overpressurized liquid (Moroni and Bonin- 



segm 



2004) 



The second is the presence in an ordinary 
helium crystal of highly mobile, extended defects. Thus, 
any manifestations of superfluid behavior of a crystal of 
helium may be inherently difficult to disentangle from 
the response of persistent pockets of liquid. 

Initial attempts to detect the possible supersolid char- 



acter of solid 4 He yielded negative outcomes (Meisel 
1992). However, some of those experiments were based 
on the misconception that a supersolid would act exactly 
like a superfluid, e.g., eliminate any pressure differences 
between two cells connected by a narrow channel ( "super- 
leak"), or allow motion of macroscopic external objects 
through the crystal. As discussed in Sec. |I.B| a true 
supersolid will only allow dissipation-free flow of its own 
constituent particles; with respect to any other object, 
its behavior is no different from that of a regular solid. 
For example, a steel ball will not sink under gravity to 
the bottom of a vessel filled with supersolid; nor will a 
pressure gradient disappear, either established within the 
supersolid sample, or between two connected chambers, 
as long as the lattice structure itself remains kinemat- 
ically stable. Of course, pressure differences will result 
in the redistribution of the superfluid component, so as 
to ensure constancy of chemical potential throughout the 
system, for a given pressure field -P(r). However, if the 
superfluid fraction is small such a redistribution will only 
result in tiny corrections to the original pressure gradi- 
ents, VP(r) — > V?(r)(l — ps/p)- This is, in essence, 
an analog of the fountain effect, but with an additional 
ingredient that the normal component is immobile even 
in the bulk (excluding plastic flow from the picture). 

At the time of this writing, there is no coherent un- 
derstanding of the observed period shift in torsion oscil- 
lator experiments. There exist a number of experimental 
observations which render problematic a straightforward 
interpretation of Kim and Chan's findings in terms of 
supersolid 4 He. For example, it is found that the period 
shift occurs in concomitance with a stiffening of the crys- 
tal pay and Beamishl 120071 iDay ei al.\ 120091); that it is 



sensitive to the resonant frequency of the apparatus; that 
it lacks the characteristic critical behavior at the alleged 
superfluid transition temperature; that it is present in 
TO experiments in which blockages are inserted, which 
should radically affect the superfluid response; that it is 
affected by the presence of small concentrations of 3 He 
impurities, as well as by DC rotation; and that it can be 
even incorrectly identified, due to the delicacy of high- 
temperature "background" subtraction (Reppy 2010). 
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FIG. 3 One-body density matrix n(r), computed by Monte 
Carlo simulation, of solid (hep) 4 He near (density p — 0.0292 
A -3 ) and away from (density p = 0.0359 A -3 ) the melting 
curve at a low temperature T = 0.2 K. Results obtained at 
lower temperature (as low as 0.1 K) are not distinguishable 
from those shown here, within the statistical uncertaities of 
the calculation. The exponential decay of n(r) signals no off- 
diagonal long-range order (i.e., Bose- Einstein Condensation) 
in the system. Reprinted from Ref. ( |Boninsegni et a/T| [2006c). 



to be reliable theoretical results for hep crystals of 4 Hc, 
and on the possible scenario of interconnected superfluid 
network in structurally disordered crystals, as a plausible 
explanation for at least some of the present experimen- 
tal evidence. This scenario is also supported by recent 
mass transport experiments by Ray and Hallock (Ray 
and Hallockl 120081). 



We shall thus focus in this paper, on what we believe 



A. Perfect hep crystals of 4 He 

Fig. [3] shows the one-particle density matrix n(r) of an 
ideal hep 4 He crystal at low temperature T = 0.2 K, com- 
puted by Monte Carlo simulations at two densities, i.e., 
near the melting curve and at relatively high pressure. By 
"ideal", we mean a commensurate crystal, i.e., one with 
exactly two atoms per unit cell; no point defects such 
as vacancies or interstitials, nor extended defects such as 
dislocations, disclinations, or grain boundaries, are in- 
cluded. These simulations are based on a microscopic 
model of helium (Eq. 12 1, making use of the accepted 
Aziz pair potential ( |Aziz et al. 1979 ) . 

The results for n(r) shown above are independent of 
temperature (below 2 K, down to at least 0.1 K), as well 
as on the number of particles in the crystal, at least for 
distances r < L/2, where L is the linear system size 
( |Boninsegni et ai\ 2006c). The observed temperature- 
independent exponential decay constitutes robust evi- 
dence that the system is insulating in the ground state 
(i.e., non-superfluid) . Such an exponential decay has 
been observed in every single calculation of the one- 
particle density matrix in commensurate solid 4 He, based 
on unbiased Quantum Monte Carlo methods. This in- 
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eludes not only finite-temperature calculations (Clark 



and Ceperley 2006 Rota and Boronat 2011), but also 



but rule it out. Moreover, effective interactions among 
vacancies 9 are found to be strongly attractive (Bonin- 



those based on ground state projection ( Galli and Reatto 
2008). It is important to note that considering a com- 
mensurate number of particles, cannot have any effect 
on the properties of the supersolid phase, in which the 
number of particles is a continuous function of the chem- 
ical potential. Insulating crystals, on the contrary, can 
change particle number in response to small changes in 
fx, at constant volume, only by adding/removing com- 
plete atomic layers. In other words they (i) always stay 
commensurate, and (ii) have zero isochoric compressibil- 
ity since the process of adding/removing atomic layers is 
kinematically frozen in an ordered state. 

A comment is in order at this point, since in this Col- 
loquium we only discuss numerical evidence such as the 
one displayed in Fig. [3] In our opinion, these results 
qualify as "proof (though not in a strict mathemati- 
cal sense) that a perfect 4 He crystal does not possess 
off-diagonal long-range order. Most numerical studies 
have well-known limitations, chiefly the finite size of the 
system that can be simulated using modern computing 
facilities. There exist, however, well-defined ways (i.e., 
finite-size scaling analysis) to deal with such limitations, 
and essentially remove them by assessing their effects 
on the physical results obtained, in the thermodynamic 
limit. It seems objectively difficult to dismiss outright 
the results of Quantum Monte Carlo simulations, in light 
of the quantitative accuracy with which they have been 
shown to reproduce many measurable properties of he- 



lium and other quantum solids (see, for instance, Ceper- 



ley (1995)). Though the validity of the prediction of ab- 



sence of ODL RO based o n simulations has been called 
into question ( [Anderson] |2009[ ), it seems fair to state 
that, given the difficulty of obtaining a formal analyt- 
ical proof, numerical results such as those quoted here 
presently constitute the strongest foundation, on which 
as unbiased as possible a statement can be made, with- 
out making any a priori assumption on the physics of the 
system. 

The exponential decay of the density matrix is consis- 
tent with finite energy gaps for creating vacancies and 
interstitials. These are readily inferred from the expo- 
nential decay of the single-particle imaginary-time (Mat- 
subara) Green function, 



G(p,r) = -(T[i> (p,r)^t(p,o)]> 



(14) 



where T is the imaginary time ordering operator and 
^(r, t) is a time-dependent field (see, for instance, Fet- 



ter and Walecka| ( |1971[ )). The above quantity was also 
computed by Monte Carlo, giving at the melting density 
Ay Ri 13 K and Aj m 23 K for vacancy and interstitial 



formation gaps, respectively ( |Boninsegni et al 2006a). 
The insulating gap A/ + Ay was found to increase with 
pressure, as expected. 

The Andreev-Lifshitz supersolid scenario is based on 
gapless vacancies or interstitials, with repulsive interac- 
tions. The large values of Ay and A/, by themselves, all 



segni et al. 2006a ), causing the formation of bound states 
of two, three, etc. vacancies. Attractive interactions 
among vacancies will lead to their collapse, the resulting 
final state depending on vacancy concentration and/or 
system density. For initial densities below freezing, i.e., 
p < pf w 0.0259 A~ 3 , the final state will be all liq- 
uid; for initial density between freezing and melting, i.e., 
Pf < p < p m ~ 0.0287 A -3 , one will observe phase sepa- 
ration into liquid and solid domains; finally, for densities 
p > Pm, the final state will contain dislocation loops pro- 
duced by coalescing vacancies. Even if out-of-equilibrium 
vacancies are introduced into the crystal at some small 
finite density, they will not be able to form a long-lived 
metastable condensate, similar to that which exists in 
cold atomic gases. The crucial difference between the two 
systems, lies in the presence of the solid matrix in helium, 
playing the role of "third body" required to satisfy energy 
and momentum conservation, in recombination processes 
leading to clustering of vacancies. 

The above results constitute strong evidence that a 
perfect hep crystal of 4 He is an insulating solid. By "per- 
fect" , we mean a single crystal, without extended struc- 
tural defects such as dislocations, disclinations, and grain 
boundaries. In turn, this suggests that, if a, hep crystal 
does have the ability of supporting dissipation-less flow 
of its own atoms, this must necessarily take place along 
extended structural defects. 



B. Superfluidity of crystalline defects 

Consider a many-body system in which the motion of 
particles is confined to one or two spatial dimensions. For 
example, one could imagine an assembly of cold atoms 
with tight optical confinement in one or two spatial di- 
rections, or a helium film adsorbed on a substrate, or a 
superconducting wire. In these cases, the rationale for 
ignoring the dimensions that are "blocked", is that the 
systems are either in vacuo, or embedded in insulating 
materials made of different atoms or molecules. When it 
comes to extended, structural defects in a crystal, how- 
ever, which can be regarded as many-body systems in 
reduced dimensions, superfluid properties are insepara- 
ble from those of the insulating bulk system. This is be- 
cause the same particles which form an insulating state 
in the bulk, are responsible for superfluidity in the lower- 
dimensional defect structure. 

Realistic systems are rarely perfect on the macro- 
scopic scale; domain walls form naturally in the out-of- 
equilibrium cooling process across the transition temper- 
ature. One may wonder whether transport properties 



9 The same results obtained for vacancies are found for interstitials 
too. Henceforth, we shall restrict our discussion to vacancies. 
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Finally, superfluidity along structural defects which 
form an interconnected three dimensional network may 
be considered as yet another scenario for obtaining the 
supersolid state. The corresponding state clearly breaks 
both translation and U(l) symmetries even though it is 
not in thermal equilibrium. Still, since extended defects 
such as domain walls, grain boundaries, and dislocations 
are protected by the lattice topology the lifetime of their 
structural relaxation may easily exceed all relevant ex- 
perimental time scales. 



FIG. 4 (a) Domain wall in the checkerboard solid, (b) When 
an atom is moved along the wall by one atomic distance, the 
wall position shifts in the perpendicular direction while the 
kink on the wall moves along the wall by two atomic distances. 



of extended defects in quantum solids are fundamentally 
different from the matrix in which they are embedded. 
An insulating state is typically the result of competition 
between kinetic and increasing potential energy terms. 
One can imagine cases when in some finite vicinity of the 
liquid-to-solid transition in the bulk system, even small 
changes in the arrangement of atoms change the balance 
in favor of particle derealization. For example, poten- 
tial barriers for moving atoms along the structural defect 
may be substantially lower than in the bulk. 

In order to gain some understanding of this physical 
mechanism, consider a model of lattice hard core bosons, 
with a repulsive interaction V between two particles oc- 
cupying adjacent lattice sites. Let —t be the value of the 
matrix element of the Hamiltonian corresponding to the 
hopping of a particle to a nearest-neighboring site. In 
a square lattice geometry and at half filling (i.e., at a 
density corresponding to one particle every two sites), a 
"checkerboard" lattice will form when V > V c = 2t. 
Fig. [4] shows a domain wall in a checkerboard crystal 



( |Burovski et al] |2005[ |Soyler et al] |2007[ ) . Such a defect 
automatically appears in this lattice model, if periodic 
boundary conditions are adopted and an odd number of 
columns is assumed, as in Fig. [4] When an atom in 
the configuration shown is shifted to a nearest neighbor 
position in the bulk, the potential energy goes up by 
3V. Shifting an atom by one atomic distance along the 
wall will increase potential energy only by V. Clearly, 
these considerations crucially depend on both the bulk 
and the defect structure, as well as the type of the motion 
considered. 

It is also easy to see from this example an intricate 
connection between the motion of atoms along the defect 
and the motion of the defect itself in the bulk. If we take 
the classical picture shown in Fig. Qa) and shift up by 
one atomic distance up all atoms nearest and to the left 
of the dashed line then the wall position will shift to the 
left in the transverse direction, see Fig. |4|h). As we will 
see shortly, this brings about interesting physics relating 
superfluidity, roughening, and quantum plasticity. 



C. Domain walls 

Most transitions into the solid phase are first-order in 
nature, and, in general, there is no reason why domain 
walls or grain boundaries in solids should have superfluid 
properties. It is the right combination of light atomic 
mass, weak inter-atomic forces, and degree of structural 
frustration, which allows atoms to establish long-range 
coherence at the defect core. For the domain wall shown 
in Fig. [4ja) one can argue, however, that at least for 
V — V c <C t, superflow of atoms along the domain wall 
will surely take place (Burovski et al. 2005 1. This is 



because for a system of hard-core bosons with nearest- 
neighbor interactions, the superfluid-solid transition is 
marginally first-order. That is, the critical point itself 
has the same extended SU(2)-symmetry as the isotropic 
Heisenberg model. 

As the interparticle interaction is increased beyond 
a threshold value Vw {V w = 3.57(3)t in d — 2 and 
V w = 2.683(3)i in d = 3 (|Soyler et al\ [2007]) one ob- 



serves the superfluid-to- insulator transition (SF-I) within 
the wall, with the concomitant disappearance of gapless 
sound excitations. In the limit of large V/t, the classi- 
cal picture of Ising-type domain wall shown in Fig. [4](a) 
is recovered. For two-dimensional walls in 3D solids, no 
obvious reason seems to exist for relating superfluidity to 
structural transformations of walls and grain boundaries, 
at the SF-I transition point. 

The same considerations apply to any domain wall, 
twin, or generic grain boundary in solids, with an impor- 
tant caveat — most solids form through a strong first-order 
transition, and the resulting crystalline structures are 
typically very efficient at localizing atoms. For example, 
even in the highly quantal helium crystals at the melting 
curve, the smallest exchange cycles are suppressed rela- 
tive to the liquid phase, by several orders of magnitude 
(despite the large zero-point motion of individual atoms) . 
As a figure of merit, one may mention the energy scale 
for the tunneling ma trix element of 3 He impurity in 4 H e 



solid, J34 ~ 10 K (Guyer et al. 1971 Richards et al. 



1975). It is more than five orders of magnitude smaller 
than the Debye frequency characterizing large zero-point 
motion of atoms, leading to an effective mass of a 3 He 
impurity which is about four orders of magnitude heavier 
than a helium atom in vacuum. In other words, struc- 
tural order emerging from first order transitions in most 
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solids is so strong, that grain boundaries and other ex- 
tended defects are likely to end up on the insulating side. 
Helium is a unique element in this regard, since it is the 
only one remaining in the liquid phase at saturated vapor 
pressure. Though experimental evidence for supe rfluidity 
of grain boundaries in 4 He is rather incomplete (Balibar 



(Shevchenko 19871 that edge dislocations in 4 He may be 



and Caupin 2008 1, Quantum Monte Carlo simulations do 



indicate that this possibility is real (Pollet et al. 2007 1 



One-dimensional walls in 2D solids allow for an in- 
teresting connection between superfluidity and roughen- 
ing, at least in some cases. Naively, one might think of 
these two phenomena as being completely unrelated, as 
they pertain to different types of motion (i.e., along and 
perpendicular to the wall direction). However, in some 
cases they may have exactly the same microscopic origin 
when neither is possible without the other. As shown in 
Fig. |4|b), the simplest, and energetically least expensive 
deformation which results in the displacement of the wall 
in the x-direction, consists of creating a "kink and anti- 
kink" pair, by shifting atoms along the wall. Clearly, 
from such a deformation mass transfer along the wall 
inevitably ensues. This suggests the intriguing possibil- 
ity that superfluidity may occur in concomitance with a 
roughening transition. By "rough" we imply a line with 
large mean square shape fluctuations (Ax) 2 (actually di- 
verging in the thermodynamic limit). More quantita- 
tively, if the instantaneous wall coordinate is described 
by x(y, t) then 



(Ax) = (£ / dy [x(y,t)-xY) , 
Jo 

r L v 

with x — L^ 1 / dy x(y,t) , (15) 
Jo 

where the average is taken over the statistical ensemble. 
The line is in the so-called "smooth" state, wherein (Ax) 2 
is system size independent. 

Interestingly enough, roughening helps to stabilize the 
superfluid phase. One might think, by looking at the 
classical picture of a straight wall shown in Fig. |4][a), 
that checkerboard order in the bulk de facto doubles the 
size of the unit cell, and therefore the effective occupation 
number in the wall is integer. Notice, however, that the 
bulk potential is shifted by half a period when the wall is 
displaced by one lattice spacing in the x-direction, mean- 
ing that for rough walls the average potential still retains 
the translation symmetry of the original lattice. Thus, 
the effective filling factor remains half-integer. This sig- 
nificantly increases the parameter range within which su- 
perfluidity survives. Numerically, one indeed finds that 



superfluid and roughening transitions coincide (Soyler 
Tied] 12007). 



superfluid, but in the absence of first-principle simula- 
tions it is impossible to predict which edge dislocations, if 
any, have that property. It is found by Monte Carlo sim- 
ulations, however, that most edge dislocations arc in fact 
insulating, their core splitting into two partials. A split- 
core configuration substantially reduces frustration in the 
atomic arrangement, and locks atoms in an insulating 
state (Pollet et al. 2008). Instead, the most basic screw 



dislocation, oriented along the c-axis of the hep structure, 
is found by simulations to support strong superfluid flow 
along the dislocation line (Boninsegni et aZ.l 120071). At 



the melting point, the effective one-dimensional super- 
fluid density is about p]f s» 1 A -1 corresponding to a 
"superfluid tube" at the core of diameter w 6 A. The 
corresponding Luttinger liquid parameter, K ss 4.9(5), is 
large enough to ensure that superfluidity of screw disloca- 
tions be stable against weak disorder or a commensurate 
potential. 

Since screw dislocations are very common in crystals, 
and often facilitate crystal growth, it is not unreason- 
able to expect that a network of superfluid "pipes" pen- 
etrating the insulating bulk may exist in real samples of 
solid 4 He, capable of conducting mass currents of 4 He 
atoms, without measurable dissipation. An obvious is- 
sue that any theoretical scenario based on a network of 
dislocations must address, is whether different types of 
dislocations with comparable superfluid properties may 
exist, in order for a fully connected three-dimensional 
network to be established. As it turns out, edge disloca- 
tions with Burgers vector along the hep c-axis and the 
dislocation core along the x-axis in the basal plane, also 
display in simulations a finite, albeit weak, superfluid re- 
sponse (Soyler et al. 2009[ ). In this case, the core is split 



into two partials separated by extremely large distance 
(possibly as large as 150-200 A) and the superfluid signal 
is "marginal", i.e., the corresponding Luttinger parame- 
ter is close to the critical value K = 2 (with a large un- 
certainty) for which a commensurate potential becomes a 
relevant perturbation. The large separation between the 
two half-cores is explained by the very small (practically 
unmeasurable) surface energy of the stacking fault defect 
in the hep structure, leading to a very weak linear poten- 
tial confining partials. The unusual split-core structure 
may have interesting consequences for how edge disloca- 
tions move, anneal, and are pinned in quantum crystals; 
these issues still need to be fully explored, at the time of 
this writing. 

Superfluidity in the cores of edge dislocations has im- 
portant consequences in terms of how solids react to 
chemical potential differences, leading to a novel phys- 
ical behavior, to which we refer as "quantum metal- 
lurgy" (though helium is obviously not a metal) 10 . At 
temperatures well below Ay, i.e., in the absence of ac- 



D. Screw and edge dislocations in He 



There are infinitely many distinct dislocation lines in 
the solid structure. It was conjectured by S. Shevchenko 



We attribute this term to A. Dorsey 
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tivated vacancies, non-superfluid edge dislocations can- 
not move in the climb direction perpendicular to their 
axis and Burgers vector. For, such a climb would corre- 
spond to the growth of an atomic layer which, in turn, re- 
quires mass transfer. Since activated vacancies is the only 
mechanism of atomic mass transport in insulating solids, 
one concludes that standard edge dislocations ought not 
climb. This, in particular, means that insulating crystals 
are isochorically incompressible, \ = {dp/dp)v ~ 0; in 
other words, the density p of the crystal does not change 
atom by atom 11 in response to infinitesimal, quasi-static 
changes in fi. 

If cores of edge dislocations in the network are super- 
fluid, then efficient transport of atoms throughout the 
sample is possible, and superfluid edge dislocations may 
to climb (or, superclimb). Correspondingly, the crys- 
tal density ought to react continuously to small changes 
in chemical potential, and an anomalous isochoric com- 
pressibility of a material which is insulating in the ab- 
sence of structural defects, should be observable. This 
description precisely matches the observation of Ray 
and Hallock, who carried out an experiment wherein 
4 He atoms (superfluid liquid) were fed into the crystal 



through implanted Vycor rods (Ray and Hallock 2008 
2009). In such a setup, the chemical potential \i is the 



physical quantity relevant to the external perturbation 
applied to the crystal. Also, the increase in the crystal 
density was correlated with the ability of solid samples 
to conduct atoms between the Vycor rods. 



E. Superglass 

As mentioned above, extended structural defects are 
metastable excited states of the system, which are pro- 
tected by topology of the crystal. One may wonder to 
what a degree this picture may continue to hold, as the 
number of defects increases, to the point where the orig- 
inal solid structure may barely be recognizable; to take 
this one step further, one may even consider a metastable, 
glassy phase arising from the application of an instanta- 
neous compression to a superfluid. 

Such a highly disordered state may be short-lived, and 
quickly relax to a polycrystal with a clearly identifi- 
able mesoscopic solid structure (here, "quickly" is defined 
with respect to relevant experimental time scale). How- 
ever, this need not be necessarily the outcome, for a glass 
produced from a sudden compression, such as described 
above, may actually be stable on the relevant experimen- 
tal (or any meaningful) time scale 12 . 



11 In the absence of vacancies and interstitials, the density of a 
crystal can react dynamically to small changes in the chemical 
potential, 8fx, only by creating or removing crystalline layers. 
This requires nucleation times exponentially large in |<5^t| — 1 . 

12 Many amorphous solids are among the strongest in Nature, with 
respect to mechanical shear or deformation. 



According to the conventional wisdom, at low tem- 
perature glasses are insulators. For, if atoms were al- 
lowed to move around, then the system would quickly 
discover a (poly)crystalline arrangement lower in energy. 
However, the underlying assumption that any kind of 
atomic derealization will eventually lead to formation 
of polycrystalline seeds, is not justified. There is ac- 
tually no fundamental reason to assume that jamming 
of structural relaxation should be incompatible with su- 
perflow. In fact, one can imagine a system where these 
two types of atomic displacements occur on vastly dif- 
ferent time scales, in such a way that the glass may re- 
main structurally stable and superfluid, under realistic 
experimental conditions. The resulting metastable phase 
breaks translation invariance, as the density profile av- 
eraged over characteristic time scales of interest is inho- 
mogeneous; it simultaneously features superfluid order. 
As mentioned in the introduction, however, since it lacks 
density LRO order, the name "superglass" (SG) is more 
appropriate than "supersolid" , for such a phase. At the 
phenomenological level, one may view the SG as a spe- 
cial limit of Tisza's two-fluid picture, wherein the normal 
fluid component freezes into a glass, instead of forming a 
quasiparticle gas. 

The first observation of SG was reported in numerical 
Quantum Monte Carlo simulations of solid 4 He samples, 
quenched from the normal liquid state at relatively high 13 
temperature an d density p = 0.0359 A~ 3 , to low temper- 
ature (~ 0.2 K) (Boninsegni et al. 2006c). A similar out- 



come is achieved when samples are prepared by squeezing 
the low-temperature superfluid liquid into a smaller vol- 
ume. When only local updates of atomic trajectories are 
employed in simulations, thus prepared samples of about 
800 atoms remain in the SG phase with inhomogeneous 
density profile and no density LRO. Though a direct con- 
nection with experiments is hard to make with regards 
to the sample preparation protocols, it seems that broad 
metastability limits for the SG state of solid 4 He are guar- 
anteed by simulations. 

Further support to the SG concept, comes from the 
possibility of "designing" a quantum system which is 



guaranteed to have a SG as an excited state (Biroli et al. 



2008). The idea is based on two key observations: 



1. Given a many-body "Jastrow" wave function for a 
system of bosons, namely an expression of the form (the 
notation is standard) 



(J) 



n 

i<j 



exp 



1 



(16) 



u being an arbitrary function of the distance = \t— Tj | 
between two particles, one can construct a Hamiltonian 
Hj, with contact two- and three-body interactions, for 



13 The reader is reminded that the density of the solid at the T= 
melting line is p m = 0.0287 A~ 3 
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which is the exact ground state wave function. One 

may check by direct substitution that is a solution of 
the Schrodinger equation for the many-body Hamiltonian 



H, 



N n 2 

2-rrii 

i—1 



+ U(r 1 



,r N 



with 



U = -- 



2m 



i<j 



i<j,i<k 



Vu(rij)-Vu(r ik ) 



(17) 



(18) 



with energy eigenvalue equal to zero (the spectrum of 
(171 is non-negative). 

2. The square of the Jastrow wave function can be 
straightforwardly interpreted as the Boltzmann-Gibbs 
statistical weight for a classical system with pair-wise 
potential V = uT. This correspondence extends to all 



eigenstates which are common for the Hamiltonian (17) 
and the Fokker-Planck operator governing the evolution 
of the classical probability distribution within the frame- 
work of the stochastic Langevin type dynamics 



1 dt 



(19) 



where r]i(t) is the thermal Gaussian white noise char- 
acterized by the correlation function (ryf (t)7/f (t')) — 
iriiSijSapSit — t'). We omit here the derivation which 
can be found in standard statistical mechanics text, see 
Ref. ( |Biroli et al.\ [2008] ) for details. 

The Jastrow wave function can describe a supersolid, 
for the following reasons: a) its square is isomorphic to a 
classical distribution, and classical many-particle systems 
have a crystalline ordered state for specific interactions 
and density. 

b) it can be shown to correspond to a quantum many- 
particle state with a finite condensate fraction. 

Quantum-to-classical correspondence can be used to 
compute both static and dynamic correlation functions at 
zero temperature, by using known classical counterpart 
results. This observation directly links the well studied 
jamming phenomenon in classical systems, to the long- 
lived metastable amorphous density profiles in the quan- 
tum case. As far as superfluidity is concerned, it is guar- 
anteed to exist by construction. Admittedly, there are 
serious drawbacks in the theoretical construction (17), 



originating from pathological properties of the resulting 
Hamiltonian Hj. This Hamiltonian always has the same 
ground state energy eigenvalue (zero), regardless of the 
system density (which implies an infinite compressibil- 
ity), a quadratic dispersion relation for elementary exci- 
tations, and a SG only appears as an excitation above 
a supersolid ground state. Nevertheless, the above ar- 
gument furnishes a proof of principle of the existence 
of SG, and can be made more physical by adding weak 
standard pair- wise potentials to Eq. (17), with the ex- 
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FIG. 5 The minimal phase defect in the network is obtained 
by imagining a vortex ring winding around the pipe and cre- 
ating ±27r phase windings in elementary plaquettes. 



F. Shevchenko state 

At the time of this writing, the most promising theo- 
retical proposal accounting for transport phenomena ob- 
served by Ray and Hallock in solid 4 He, seems to be the 
3D network formed by interconnected one-dimensional 
superfluid channels. Let us start the discussion with a 
fictitious system of pipes running along bonds of a sim- 
ple cubic lattice, see Fig. [5j Assume that pipes are of 
diameter d=l mm, filled with low-temperature (T <C T\) 
Helium-II, and the bond length L is such that L ^> d /a 
where a is the interatomic distance, say L = 10 6 km. 
The reader might think that the thermodynamic tran- 
sition temperature to the SF phase, in this quasi one- 
dimensional setup, should be strongly suppressed, rela- 
tive to the A-point in bulk Helium-II. However, it is ob- 
vious that at, say, T < T\/2, this network will support 
frictionless, low-velocity flows, with persistent currents 
whose characteristic decay times exceed the age the Uni- 
verse! Indeed, at this temperature Helium-II already has 
strong local order (definitely on a millimeter scale) in 
the phase field, which prevents large vortex excitations 
from nucleating and proliferating in the space occupied 
by helium. 

Both expectations are correct, and not in conflict with 
each other. Simply, one is dealing here with an "extreme" 
case, wherein one kinetic time scale rapidly becomes so 
large, that in an actual experiment the ergodic hypothe- 
sis, upon which all equilibrium thermodynamic ensemble 
calculations are based, is violated. To make this state- 
ment more quantitative, consider the kinetic energy of 
a current-carrying state corresponding to a phase wind- 
ing worth 27T around one elementary plaquette, regarded 
as isolated from the rest of the system. A simple cal- 
culation based on the integration of the energy density 
e = ( Ps /2m){V<p) 2 = (p s /m)(^ 2 /8L 2 ) gives 



pectation that the overall picture may not change. 



E, 



pi 
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which is much smaller than T\ for L ^> d 2 /a (for pa- 
rameters mentioned above this energy is in the micro 
Kelvin range). In the 3D network, currents are not con- 
fined to elementary plaquettes. The minimal phase defect 
which can be imagined is a small "vortex ring" type con- 
figuration depicted in Fig. [5] Though vortexes are not 
even defined on length scales smaller than ~ L, and the 
best microscopic description is provided by phase gradi- 
ents along bonds, the topology of phase windings is still 
easiest to picture using hypothetical vortex lines. Nev- 
ertheless, Eq. (20) sets the scale for the energy cost of 
circulating currents in the network. As long as E v \ <C T 
the thermodynamic equilibrium state remains normal be- 
cause on all scales the closed-contour integrals § (V(f) ■ dr 
are taking non-zero values. 

The normal-to-superfluid transition temperature can 
be estimated from the condition E v \ ~ T, at which point 
plaquette currents become thermodynamically unfavor- 
able. This condition leads to 



Ta^oc^ (L»d 2 /«) 
La L 



(21) 



The same estimate for T c follows from the condition that 
ID phase fluctuations on the length-scale L are reduced 
to a value of order unity: by introducing an effective 
one-dimensional superfluid stiffness A^ 1D ^ ~ p$d 2 jm we 
find that ((^(0) - <p(l)) 2 ) ~ TL/A S . Once phase fluc- 
tuations along the pipe length become small the three- 
dimensional order in the network sets in. The univer- 
sality class of the transition is not altered by the quasi 
one-dimensional microscopic geometry; in particular, the 
ps(T) curve starts with an infinite derivative at T c . By 
counting the amount of superfluid liquid in the pipes we 
also find p s (T = 0) - l/L 2 ~ T 2 . 

Clearly, the thermodynamic transition discussed above 
is of purely academic interest, i.e., one does not expect to 
observe it experimentally. At temperatures below T\ , the 
notion of thermal equilibrium quickly becomes irrelevant 
because any given distribution of persistent currents is 
kinetically frozen and cannot change in response to tem- 
perature variations or slow rotation. In other words, in 
terms of observed behavior, the state is indistinguishable 
from the genuine superfluid, despite the fact that typ- 
ically one finds large circulating currents on most pla- 
quettes. This is the essence of the state introduced by 
S. Shevchenko (IShevchenkol 119871) to describe the possi- 



ble superfluid properties of a dislocation network in the 
temperature interval T c -C T -C T\, under the hypothesis 
that edge dislocations in 4 He feature superfluid cores. 

An important quantitative difference between the fic- 
titious "pipelines" example and dislocation network, is 
that in the latter the pipe diameter is of the order of the 
interatomic distance, d ~ a. The condition L 3> d 2 /a 
is modified into L 3> a, and thus is easily satisfied un- 
der realistic experimental conditions except, may be, in 
a glassy phase. Correspondingly, one may expect that 
T c -C 7a, i.e., that there should be a broad tempera- 
ture interval where the Shevchenko state may be real- 
ized. Another distinction is that narrow superfluid pipes 



at low-temperature should be described by the Luttinger 
liquid theory with relatively small values of K, i.e., their 
superfluid response is rather fragile, with stability of cir- 
culating currents being no longer protected by the ro- 
bust superfluid order on large length scales across the 
channel. Long relaxation times, tll(T), for equilibra- 
tion/redistribution of plaquette currents leading to the 
Shevchenko state within the Luttinger liquid theory are 
described by the power-law dependence (not mentioning 
back-scattering matrix elements leading to momentum 
non-conservation) 



tll(T) 



2K-1 



(22) 



where T + ~ mc 2 is the characteristic energy scale based 
on the one-dimensional sound velocity c below which the 



Luttinger liquid behavior sets in (Kashurnikov et al. 
1996). For dislocations in solid helium T* ~ T\. At 
temperatures T ~ T\, one finds relatively short relax- 
ation times ensuring normal state behavior at experimen- 
tal time scales r oxpor . However, at low temperatures, say 
for K rj 5 found for screw dislocations, it is possible to 
observe a dynamic crossover to the superfluid state when 
t ll{T) ^> r cxper . The crossover might be rather sharp 
and easily accessible experimentally for large values of 
K. Note also that typical values of the Luttinger param- 
eter for ultra-cold atomic systems in the absence of the 
optical lattice are enormous, of the order of a hundred or 
larger ! The question of dynamic crossover did not arise 
for the fictitious "helium pipelines" system because the 
corresponding Luttinger parameter is proportional to the 
channel diameter squared K ~ {d/a) 2 and the phase-slip 
dynamics is frozen right below the A-point. 



III. SUPERSOLIDITY IN ULTRACOLD ATOMIC 
SYSTEMS 

Regardless of how the current controversy over the in- 
terpretation of the present 4 He experiments is eventually 
resolved, it seems fair to state that solid helium does not 
afford a direct, simple, and clear observation of the su- 
persolid phenomenon. The question then arises of which 
other physical system may allow one to make a rela- 
tively easy, unambiguous experimental identification of 
this novel phase. 

Among all simple atomic or molecular condensed mat- 
ter systems, helium offers by far the most favorable 
combination of large quantum derealization of its con- 
stituent (Bose) particles, owing to the light mass of its 
atoms, and weakness of the interatomic potential. It is 
precisely for this reason that liquid helium escapes crys- 
tallization at low temperature, under the pressure of its 
own vapor. The closest condensed matter system that 
may enjoy similar properties, is molecular hydrogen (H2), 
also an assembly of Bose particles. Indeed, the mass of 
a H 2 molecule is one half that of a helium atom, which 
would lead one to expect even higher quantum effects. 
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However, the attractive well of the interaction between 
two hydrogen molecules is about three times deeper than 
that of two helium atoms. As a result, liquid hydrogen 
crystallizes at a relatively high (14 K) temperature, sig- 
nificantly above that at which BEC and SF are expected 
to occur; although quantum effects, including those of 
Bose statistics, are detectable in the momentum distri- 
bution of the liquid near melting ( Boninsegnil 120091) , and 



although superfluid (and even supersolid) behavior has 



been predicted for small clusters of parahydrogen (Mez 



zacapo and Boninsegni 2006 2007 2011 ), in general, the 



behavior of solid molecular hydrogen is much closer to 
that of a classical crystal than to solid helium. No exper- 
imental evidence for possible super fluid behavio r of so lid 
hydrogen has so far been reported ( |Clark et al. 



2006). 



The behavior predicted for solid 4 He, as it emerges 
from first-principle quantum simulations, is largely de- 
termined by the strong repulsive core of the interatomic 
potential at short distance (less than ~ 2 A). While 
the attractive (Van der Waals) long-range part of the 
interaction is responsible for the existence of the con- 
densed phase, it is the repulsive core (which is a result 
of Pauli exclusion principle, acting so as to prevent elec- 
tronic clouds of different atoms from overlapping) that 
determines most of the thermodynamic properties of con- 
densed helium and other quantum solids and liquids. For 
example, a very simple model of Bose hard spheres re- 
produces surprisingly accurately the phase diagram of 
condensed helium. 

Computer simulation studies of classical crystals, mak- 
ing use of the Lennard- Jones potential, have yielded ev- 
idence of the same vacancy phase separation observed 
in the quantum system ( Ma et al. 2008 1 . This suggests 



that the origin of the thermodynamic instability of a gas 
of point defects lies in the strong interaction among par- 
ticles, which quantum derealization cannot overcome. 
One is therefore led to consider systems characterized by 
a different type of pair-wise interaction, possibly with 
a "softer" core at short distance. The question is, of 
course: Where, in nature, does such an interaction arise? 
The basic features of the helium interatomic potential, 
chiefly the strong repulsion at short distance, are com- 
mon to any molecular or atomic interaction. 
One route consists of searching for the supersolid phase in 
systems in which the "elementary" constituents are com- 
posite particles, e.g., Cooper pairs in superconducting 



Josephson junction arrays (Roddick and Stroud 1995 1, 



or excitons in electron-hole bilayers Ooglekar et 



al. 



2006 ). In this case, particles dynamically form and disap- 



pear due to pair breaking and recombination effects; the 
interaction between two such objects is of the effective 
kind, i.e., induced by the medium in which these parti- 
cles are embedded, and has an important time-dependent 
component. On the other hand, if one wishes to retain 
the simple picture of "elementary" particles interacting 
via a static pair potential, dilute assemblies of spatially 
confined ultracold atoms appear to offer a viable option. 
Seventeen years after the first successful observation of 



BEC in a spatially confined assembly of Rubidium atoms, 
cooled down to temperature in the nanoKelvin range 



(Anderson et al. 1995 Pethick and Smith 



2005 ) , impres- 



sive scientific and technological advances have made the 
field of ultracold atoms the natural playground, where 
fundamental issues in condensed matter and many-body 
systems can be addressed. Dilute assemblies of cold 
atoms constitute an almost ideal many-body system o) 
virtually free from the imperfections and "background 
noise" that mask subtle physical effects in a real solid, 
and b) upon which a remarkable degree of control can 
be achieved. In particular, there exist a number of 
techniques allowing one to vary the interaction between 
atoms or molecules, rendering it essentially an adjustable 
parameter. The simplest, and so far most commonly 
utilized such technique, takes advantage of the so-called 
Feshbach resonance (see, for instance 
(2005)), whereby the strength of a (: 



Pethick and Smith 
lort-ranged) inter- 
action between two atoms or molecules can be varied, 
and even its sign reversed (i.e., turned from repulsive to 
attractive, or vice versa). 

"Fashioning" artificial interparticle potentials, not 
arising in any known condensed matter system, allows 
one to address a key theoretical question, namely which 
two-body interaction potential(s), if any, can lead to the 
occurrence of a supersolid phase. 

In the next two subsections we shall discuss two physi- 
cal systems, both realizable in the laboratory using cold 
atoms or molecules, which may provide a direct pathway 
to the stabilization and presumably straightforward ob- 
servation of a supersolid phase, owing to the particular 
interactions among elementary constituents. 



A. Rydberg Blockade and soft-core potentials 

TheRydberg blockade is a physical mechanism that was 
initially introduced as a way to manipulate quantum in- 
formation stored in collective states of mesoscopic en- 
sembles (Lukin et al. 2001), has been recently proposed 



as a way to engineer a novel type of interaction poten- 
tial between cold atoms. Specifically, the modified in- 
teraction "flattens off" , and remains essentially constant 
below some characteristic "cut-off" distance a ( |Henkel| 



et al. 2010). In particular, it is v(r — > 0) — V, finite 



and not much larger than the characteristic energies at 
play 14 . 

What is the connection between such a potential and 
the supersolid phase ? 

Within the classical complex field description based on 
the Pitaevskii-Gross equation, the answer was first pro- 



14 For comparison, one may note that the most accurate model 
potential between two helium atoms l |Aziz et al.\ |1979| is also 
finite in the r — > limit, but its value, of the order of 10 6 K, 
greatly exceeds any energy scale relevant to ordinary condensed 
matter physics. 
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(a) 



(b) 



FIG. 6 One-dimensional crystal of particles interacting via a 
soft-sphere potential, equal to some energy V > if the dis- 
tance between two particles is less than their effective diame- 
ter a, and zero if it is greater. As the density equals 1/a, con- 
figuration (b) becomes energetically advantageous over (a), as 
its cost is V/2 per particle, as opposed to V as in (a). 



vided by E. Gross more than half a century ago 



1957 ), and more recently quantified by Josserand et al. 



Gross 



( Josserand et al. 2007 ) who observed that at sufficiently 



high density, the system will break translation invariancc 
and develop solid LRO. As in any other classical field 
model, the ground state is necessarily superfluid, since 
the lowest energy configuration corresponds to perfectly 
ordered phases. Note that the notion of a "particle" is 
completely lost within the classical field description, and 
thus this study corresponds to the limit of infinite number 
of particles (and vacancies) per unit cell. Such a treat- 
ment is tantamount to regarding a supersolid as a su- 
perfluid with a density modulation. Similar results were 



obtained in Ref. (Hcnkel et al. 2010) 



Consider for definiteness a simple model in which par- 
ticles behave as soft spheres, namely interact via a "box" 
potential, equal to some energy V > if two particles 
are at a distance r from each other less than their ef- 
fective diameter a, and zero otherwise. These above 
considerations suggest that, if V is small, and the num- 
ber of particles per unit cell is large, a supersolid phase 
will be stable if the dimensionless coupling parameter 
g = (Vma 2 /Ti 2 )(pa?) is greater than unity (not too large, 
or the system will behave essentially as a classical crys- 
tal). At the microscopic level, one can intuitively under- 
stand why, as a consequence of the Constance of the po- 
tential at short distance, at sufficiently high density the 
system will find it energetically favorable to form crystals 
with relatively large numbers of particles per unit cell. 
As a simple illustration, consider the one-dimensional 
crystal shown in Fig. [6]ja). Assume again the simple 
soft-sphere model described above. Clearly, at a density 
1/a each particle can be made to touch just two oth- 
ers, as shown in Fig. [6ja), at a total energy cost of V 
per particle. The system can lower its potential energy, 
by slightly displacing every other particle (i.e., moving 
pairs of particles closer together) and creating a crystal 
of lattice constant 2a, with two particles per unit cell, 
as shown in Fig. |6jb). As the density is increased, the 
number of particles per units cell also increases. 
The outstandin g questions, v i s-a-vi s the predictions 
made in Refs. (Henkel et al., 2010[ Josserand et al. 
[2007] ), are a) how large a number of particles per unit 
cell one might require in practice, and 6) what would be 
the physical nature of such a supersolid phase. 



% ♦ 1§ ; -if 




FIG. 7 Snapshots from computer simulations of a system 
of bosons interacting via a power-law potential flattening off 
at short distances (see Ref. ( Cinti et al. \ |2010[ ) for details) . 
Snapshots are taken at four different temperatures, decreas- 
ing from (a) to (d); each snapshot is taken at a temperature 
an order of magnitude lower than the immediately previous 
one. Reprinted from Ref. (Cinti et al. 20101. 



In order to have a fully quantum mechanical quanti- 



tative solution of the problem Cinti et al. (Cinti et al. 



2010 ), investigated by Quantum Monte Carlo simulations 



the low temperature phase diagram of a two-dimensional 
assembly of Bose particles, interacting via a number of 
different potentials of the general form v(r) ~ l/r™ for 
r > a, but either equal to a constant or slowly vary- 
ing for r < a, in any case smoothly approaching a finite 
value for r — > 0. In those studies, n was taken equal to 
3 and 6. As it turns out, the main physical results are 
independent on the behavior of the potential for r > a. 
Indeed, recent work by Saccani et al. shows that the 
simple model of "soft spheres" described above, in two 
dim ensions, features the same basic physics outlined be- 
low (Saccani et al. 2011). 



A typical result is shown in Fig. [7] Displayed are snap- 
shots of Feynman's trajectories in imaginary time from a 
computer simulation. In the particular case shown, the 
average intcrparticlc distance is approximately a/2, and 
the values of the temperature (which decreases from (a) 
to (d)) are such that each is an order of magnitude higher 
(lower) than the next (previous) one (see Ref. ( |Cinti| 
et al.\ |2010 ) for details). As T is decreased, particles 
bunch into mesoscopic droplets, in turn forming a regular 
(triangular) crystal. We henceforth refer to this phase as 
the droplet- crystal phase. The formation of such droplets 
is a purely classical effect; indeed, simple (but remark- 
ably accurate) estimates of the average number N4 of 
particles per droplet, can be obtained based on purely 
classical potential energy minimization, for the various 



18 



potentials considered. 

In the T — > limit, long exchanges of identical par- 
ticles can take place, as a result of particles tunneling 
from one droplet to an adjacent one. In turn, long ex- 
changes of particles can result in a finite superfluid re- 
sponse throughout the whole system, and indeed a bulk 
superfluid signal is observed at low temperature, roughly 
in correspondence of an inter-particle distance r s of the 
order of (slightly less than) the potential cutoff radius. 
The number Nd of particles per droplet for which a super- 
solid phase is observed, is variable, but can be as low as ~ 
4, for specific choices of r s and a. Because superfluidity 
arises in concomitance with the droplet-crystal structure, 
the denomination supersolid seems indisputable in this 
case. Supersolid behavior in this system originates from 
tunneling of particles between droplets which are them- 
selves individually superfluid, as simulation results show. 
This is reminiscent of the phase-locking mechanism in a 
(self-assembled) array of Josephson junctions. 



B. Dipolar systems 

The supersolid phase described in the previous section 
is a direct consequence of the "flatness" of the poten- 
tial at short distance, while the long-range behavior is 
largely irrelevant. On the other hand, supersolid behav- 
ior can also be underlain by long-ranged interactions. We 
consider here the case of atoms or molecules possessing a 
finite electric dipole moment. These particles can be con- 
fined to quasi-2D, by means of an external harmonic po- 
tential in the direction perpendicular to the motion (the 
so-called "pancake" geometry) . Upon aligning all dipoles 
in the direction perpendicular to the plane, by means of a 
strong electric field, one can study a system of Bosons in- 



teracting via a purely repulsive potential ( Buchler et al. 
2007[) of the form 1/r 3 . 



What renders such a system particularly intriguing, is 
the existence of an exact theoretical result ( |Spivak and 
Kivelson[ 2004), excluding the occurrence of first-order 
phase transitions involving a density change Sn (such as 
liquid-solid), in 2D systems with such an interaction. The 
reasoning goes as follows: a first-order phase transition is 
characterized by the coexistence of two phases of different 
density, separated by a macroscopic interface. However, 
a straightforward calculation shows that the energy of 
such an interface contains a negative term, which diverges 
logarithmically in the thermodynamic limit. 

An immediate consequence of the above result is that, 
on approaching the transition from the low-density (liq- 
uid) phase, the system will lower its free energy by em- 
bedding sufficiently large solid domains (i.e., macroscopic 
"bubbles") inside the liquid. At T — two effects are ex- 
pected to occur, namely a) the transition of the liquid to 
a superfluid, and 6) the crystallization of solid bubbles 
into a lattice superstructure, resulting in a global super- 
solid phase (in fact, a whole set of different such phases 
(|Spivak and Kivelson|[2004|). 



An intrinsic subtlety of this scenario is the com- 
petition between the positive contribution to the sur- 
face tension, a + , originating from short-range physics, 
and the negative scale-dependent contribution er_ (R) cx 
— (Sn) 2 h^Rp 1 / 2 ), where R is the droplet size. For small 
density differences (and quantum Monte Carlo seem to 
indicate that this is indeed the case) the relevant length 
scale for having negative surface tension er(i?) = u + +cr_ 
might be astronomically large and outside the reach of 
realistic experimental/numerical setups. 



IV. CONCLUSION 

The investigation of the possible supersolid phase of 
matter in helium is ongoing, and at the time of this writ- 
ing no single theoretical framework has emerged as the 
accepted interpretation of the puzzling and controversial 
phenomenology. In this Colloquium, we have presented 
one possible scenario, which is based on the superfluid 
properties of extended defects, chiefly dislocations. This 
scenario is based on first principle microscopic calcula- 
tions, whose main quality, in our view, is the lack of any a 
priori assumption on the physical behavior of the system. 
It seems consistent with at least an important part of the 
phenomenology, while attempting to establish a connec- 
tion between it, and microscopic mechanisms. Whether 
it will stand the test of time hinges on its ability to of- 
fer quantitatively more accurate microscopic predictions, 
e.g., on the supersolid transition temperature, if it is ul- 
timately established that what has been observed in the 
laboratory, is indeed a signature of superfluidity of the 
helium crystal. It should be mentioned that, since its 
introduction, there has been further, independent theo- 
retical work on the Shevchenko state (or, very similar the- 
oretical pictures) by other authors who adopted a differ- 



2008) 



ent, macroscopic approach (Goswami et al. 2011 Toner 



We have also discussed the possibility of realizing 
the supersolid phase in a different context, namely cold 
atoms. This is motivated by our belief that such a physi- 
cal setting may well afford the unambiguous experimental 
observation of this intriguing phase, more directly (and 
with less controversy) than in solid helium. 
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